Nonlinear atom optics and bright gap soliton generation in finite optical lattices 
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We theoretically investigate the transmission dynamics of coherent matter wave pulses across 
finite optical lattices in both the linear and the nonlinear regimes. The shape and the intensity 
of the transmitted pulse are found to strongly depend on the parameters of the incident pulse, in 
particular its velocity and density: a clear physical picture for the main features observed in the 
numerical simulations is given in terms of the atomic band dispersion in the periodic potential of the 
optical lattice. Signatures of nonlinear effects due the atom-atom interaction are discussed in detail, 
such as atom optical limiting and atom optical bistability. For positive scattering lengths, matter 
waves propagating close to the top of the valence band are shown to be subject to modulational 
instability. A new scheme for the experimental generation of narrow bright gap solitons from a wide 
Bose-Einstein condensate is proposed: the modulational instability is seeded in a controlled way 
starting from the strongly modulated density profile of a standing matter wave and the solitonic 
nature of the generated pulses is checked from their shape and their collisional properties. 
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In recent years, a great interest has been devoted to 
theoretical as well as experimental studies of the propa- 
gation of matter waves in the periodic potential of opti- 
cal lattices. The first experiments were carried out us- 
ing ultracold atomic samples Q); later on, the realiza- 
tion of atomic Bose-Einstein condensates (BECs) 0j and 
their coherent loading into optical lattices ||] have opened 
the possibility of investigating features which follow from 
the coherent nature of the Bose-condensed atomic sam- 
pie §§. 

At the same time, the propagation of light waves in lin- 
ear and nonlinear periodic dielectric structures has been 
a very active field of research: global photonic band gaps 
(PBG) have been observed || and one-dimensional non- 
linear periodic systems such as nonlinear Bragg fibers 
are actually under intense investigation given the wealth 
of different phenomena including optical bistability, mod- 
ulational instability and solitonic propagation that can 
be observed |,|, 0. 

Given the very close analogy between the behavior of 
coherent matter waves and nonlinear optics, we expect 
that the concepts currently used to study the physics of 
nonlinear Bragg fibers can be fruitfully extended to the 
physics of coherent matter waves in optical lattices: the 
optical potential of the lattice plays in fact the role of the 
periodic refractive index, the atom-atom interactions are 
the atom-optical analog of a Kerr-like nonlinear refrac- 
tive index and the Gross-Pitaevskii equation of mean- 
field theory corresponds to Maxwell's equation with a 
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nonlinear polarization term 

In the present paper, we shall report a theoretical in- 
vestigation of the transmission dynamics of coherent mat- 
ter pulses which incide onto a finite optical lattice with 
a velocity of the order of the Bragg velocity. In this ve- 
locity range, Bragg reflection processes are most effective 
and the atomic dispersion in the lattice is completely dif- 
ferent from the free-space one. Depending on the value 
of the density, spatial size and velocity of the incident 
atomic cloud, as well as on the depth and length of the 
lattice, a number of different behaviors are predicted by 
numerical simulations; here, we shall focus our attention 
on the shape of the transmitted pulse just after it has 
crossed the lattice as well as while it is still propagat- 
ing in the lattice. In particular, we shall discuss a new 
mechanism which can be used to generate narrow bright 
atomic gap solitons propagating along the lattice start- 
ing from a wide incident Bose-Einstein condensate. For 
more details on the continuous-wave transmission and 
reflection spectra at linear regime, the reader can refer 
to []l2[ |b| [14| [l5| |l6) ; some aspects of the linear pulse 
dynamics are discussed in 15 1. 

The geometry considered in the present paper as well 
as in Jl2], [l3| [lj, [H| |l6| is significantly different from the 
one usually considered in recent experimental works Q 
of BEC dynamics in optical lattices in which an optical 
lattice is switched on and superimposed to a stationary 
condensate; the dynamics of the condensate inside the 
lattice is then studied in response to some external force 
such as gravity, an acceleration of the lattice, or a spatial 
translation of the magnetic potential. 

The present paper is organized as follows: sec. I de- 
scribes the physical system under examination. In sec.|| 
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we review some basic concepts on the dispersion of mat- 
ter waves in the periodic potential of an infinite optical 
lattice and present some simple analytical calculations 
which accurately reproduce the dispersion of the atomic 
bands in the neighborhood of the first forbidden gap. 
The linear regime propagation of coherent matter wave 
pulses across finite lattices is the subject of sec. [II. The 
intensity as well as the shape of the transmitted pulse 
are found to strongly depend on the properties of the 
weak incident pulse, in particular its velocity and spatial 
size; a simple interpretation of the observed phenomena 
in terms of allowed bands and forbidden gaps is provided. 
In sec. IV, we discuss the effect of atom-atom interactions 



on the propagation of the pulse in the different cases: 
simple explanations for the observed behavior are put 
forward in terms of familiar concepts of nonlinear optics 
such as optical limiting, optical bistability, or modula- 
tional instability; in particular, modulational instability 
is shown to occur for positive scattering length atoms 
when the matter wave is propagating at the top of the 
valence band; at this point of the dispersion, the effec- 
tive mass is in fact negative. In sec.[y|, we propose a new 
scheme to control the modulational instability of valence 
band atoms in order to obtain a narrow bright gap soliton 
from a wide condensate: as an initial seed for the insta- 
bility, the standing matter wave pattern created by the 
interference of the incident and reflected waves is used. 
The solitonic nature of the generated pulses is verified 
by looking at their dynamical and collisional properties 
as well as by comparing the pulse shape with the an- 
alytical predictions o f th e envelope-function approxima- 
tion discussed in sec. VI, Conclusions are finally drawn 



in sec. VII 



I. THE PHYSICAL SYSTEM 

We consider a Bose-condensed atomic cloud in a quasi- 
1D geometry in which the transverse motion is frozen by 
the confining potential of an optical or magnetic atomic 
waveguide | p7[ . Gravity is made immaterial either by 
placing the waveguide axis along the horizontal plane or 
by counterbalancing the gravitational field with a suit- 
able magnetic field gradient. 

A periodic potential is created along the waveguide 
axis by means of a pair of far-off resonance laser beams 
of frequency u)l and wave vector — lul/c cross- 
ing the waveguide at an angle 9 as in fig|l] @: denot- 
ing with £Il(z) = \d ■ E(z)\/H the (slowly varying) sin- 
gle beam Rabi frequency and with u> at the atomic tran- 
sition frequency, the optical potential experienced by 
the atoms is given by V op t(z) = Vq(z) cos 2 /cbi-z, with 
Vq(z) = MIl(z) 2 /(u>l — oj at ) and fce r = /c^cos^. For 
a red- or blue-detuned laser field, the optical potential 
is respectively attractive or repulsive; the lattice period 
£bi = Tr/kBr can be tuned by varying the angle 9. The 
longitudinal envelope Vq{z) of the lattice potential is de- 
termined by the profile of the laser beam waist and is 
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FIG. 1: Upper panel: schematic plot of the experimental 
set-up under consideration. Lower panel: spatial dependence 
of the local band edge energies. 



assumed to smoothly vary on a length scale significantly 
longer than the lattice period ^bi-; unless differently spec- 
ified, Vq(z) will be taken as a Gaussian 



V (z) = V e- z2 / 2v 



(1) 



of height Vb and spatial length wi . 

If both the kinetic and the interaction energy are much 
smaller than the transverse level spacing of the waveg- 
uide, this latter can be considered as being a single-mode 
one and the condensate wave function can be written in 
the factorized form 



V>(x) = 1p(z)lpx(x.±) 



(2) 



where i/j±(x±) is the ground-state eigenfunction of 
the transverse confining potential with the appropriate 
Jd 2 x± \t/j±(x±)\ 2 = 1 normalization. Within this ap- 
proximation, the dynamics of the condensed atomic cloud 
can be described by a one-dimensional Gross-Pitaevskii 
equation 
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where mo is the atomic mass and the renormalized ID 
effective interaction <?id is written in terms of the usual 
3D scattering length a as fl8| 
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II. 



ALLOWED BANDS AND FORBIDDEN 
GAPS 



As it happens to electrons in crystalline solids Jl9] and 
light in periodic dielectric systems such as photonic band 
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gap crystals M, the atomic dispersion in an infinite pe- 
riodic potential is characterized at linear regime (i.e. in 
the non-interacting case) by allowed bands and forbidden 
gaps. 

When the depth Vb of the lattice potential is weak 
with respect to the Bragg energy Hlubi — F^k^/lmo the 
lowest part of the band structure can be accurately de- 
scribed within a nearly-free atom approximation |lq| in 
which only two coupled modes [|7[ are taken into account 
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In this restricted (f , b) basis, the linear regime Hamilto- 
nian has the following form 
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and the eigenenergies hu!± are equal to 

flU±(k) = htUBr + y + 
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where we have set Afc = k — fcB r ; the positive sign holds 
for the upper, conduction band and the minus sign for 
the lower, valence band. No states are present between 
^Br + Vo/4 and wsr + 3Vb/4: this is the lowest forbidden 
gap in which the matter waves can not propagate through 
the lattice. 

The group velocity is given by 
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close to band edge (Hvbt Afc <c Vb), the group velocity 
Vg is much reduced with respect to the Bragg velocity 
«Br = SJcBr/mo; further away (Hvbi Afc ^ Vb), it recovers 
the free-space value hk/m. 

The effective mass of the atom is given by 
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At band edge (Afc = 0), m c s is much smaller in absolute 
value than the free-space atomic mass m ; as usual, it is 
positive at the conduction band edge, while it is negative 
at the valence band edge. Further away from the band 
edge, TO c ff coincides with the positive free-space value 

TO 



Around the band edge, the weights of the forward and 
backward traveling waves are comparable |af| ~ |ab| for 
both valence and conduction bands: the density profile of 
the Bloch eigenfunction thus has a standing wave shape 
with a spatial period equal to the lattice period £bt- Fur- 
ther away from the gap one component dominates the 
other: the Bloch eigenfunction has then a running wave 
character and the density profile is uniform over the unit 
cell of the lattice. 



III. PROPAGATION IN THE LINEAR REGIME 

We now consider a coherent matter pulse (e.g. ex- 
tracted from a Bose-Einstein condensate) which is mov- 
ing along the waveguide with a uniform velocity vq close 
to the Bragg velocity «Br- Initially the atomic pulse is 
far outside the lattice. The initial density distribution of 
atoms in the cloud is taken as having a Gaussian shape 
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the carrier momentum is hko — too^o and the carrier 
kinetic energy Hujq — h 2 k^/2mQ] since the wave packet 
is finite in space, its Fourier transform tp{k) has a finite 
momentum spread o~k — 1/wq around fco- 
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FIG. 2: Upper (a) panel: linear transmission spectrum across 
a repulsive (Vo/uobi = 1) optical lattice with Gaussian profile 
(wi/£bi = 5). Lower panels: linear regime time-evolution of 
Gaussian incident pulses (wo/Ibt = 20) centered at respec- 
tively u>o — u) a (b) (transmitted) and luq = tup (c) (reflected). 
The spatial extension of the lattice is indicated by the vertical 
dashed lines. 

If the density is sufficiently low and the interac- 
tions sufficiently weak, the nonlinear term in the Gross- 
Pitaevskii equation (||) can be neglected and the time- 
evolution of the matter field is accurately described by a 
linear Schrodinger equation. In this case, the superpo- 
sition principle holds and the transmission of the pulse 
can be well described in momentum space in terms of the 
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wave vector-dependent transmission amplitude t(k) 

#t(k) = t(k)# ims (k). (ii) 

If the whole of the incident wave packet is contained in 
either a transmitting or a reflecting region of the spec- 
trum (fig.||a), it will be respectively transmitted (fig.||b) 
or reflected (fig.|]c) without sensible reshaping. As it has 
been discussed in full detail in previous papers jl6) , com- 
plete transmission occurs whenever propagating states 
are available at all spatial positions for the frequencies 
of interest. Given the smooth envelope of the lattice, in- 
terface reflections do not occur and the matter wave adi- 
abatically follows the shape of the corresponding Bloch 
state; the typical sinusoidal shape of band edge Bloch 
wave functions j^, multiplying the broad pulse enve- 
lope can be clearly recognized while the wave packet is 
crossing the lattice (fig.|]b, inset). On the other hand, if 
the wave packet frequency falls inside the forbidden gap 
at some spatial positions, the matter wave is not able to 
cross that region and is then nearly completely reflected 
back (fig.||c); the inversion point is located at the begin- 
ning of the forbidden region, i.e. at the point at which 
the carrier frequency loq enters the forbidden gap. 
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FIG. 3: Lower panel (c): time-evolution of a weak and short 
(wo /for = 10) Gaussian pulse close to resonance with a quasi- 
bound mode at lj 7 ~ 1.17cje r (see fig.^a); lattice parame- 
ters as in figs.[l] and [| The spatial extension of the lattice 
(wi /£bi = 5) is indicated by the vertical dashed lines. Upper 
panels: (a) expanded view of the transmitted pulse shape for 
the pulse in (c) (solid line) and for a longer wo/£bt ~ 20 pulse 
(dashed line); (b) incident (dashed) and transmitted (solid) 
spectra for the incident pulse in (c). 

When one or more resonance peaks due to quasi-bound 
modes 1 are contained within the spectrum of the incident 



wave packet, the pulse will be partially reflected and par- 
tially transmitted (ngj^b-c); the strong frequency depen- 
dence of the transmission amplitude leads to a profound 
reshaping of the spatial profile of the pulse (fig.^a) . An 
incident pulse of linewidth much wider than the quasi- 
bound mode linewidth has in fact a time duration much 
shorter than the characteristic decay time of the mode; 
this latter can therefore be considered as being impul- 
sively excited by the incident pulse and than exponen- 
tially decaying on a much longer time scale. The expo- 
nential tail shown by the transmitted pulse profile when a 
single quasi-bound mode is excited (dashed line in figj|a) 
is a clear signature of transmission occurring via a single 
resonant quasi-bound mode j2(J. 

If several isolated modes are instead present in the fre- 
quency interval encompassed by the incident spectrum, 
the spectrum of the transmitted pulse will contain sev- 
eral peaks (fig.||c) and a complex shape will result from 
the interference of the different frequency components. 
For example, if the incident pulse has a duration com- 
parable to the dephasing time of a pair of neighboring 
modes (i.e. the inverse of their frequency difference), 
both of them will be impulsively excited and the profile of 
the transmitted pulse will show oscillations following the 
time evolution of the relative phase (solid line in fig.[3]a) ; 
these oscillations can be interpreted as quantum beatings 
between the two coherently excited quasi-bound modes 
coupled to the same continuum of propagating modes 
outside the lattice 2 



IV. NONLINEAR REGIME AND 
MODULATIONAL INSTABILITY 



For higher values of the atomic density and the cou- 
pling constant <7id, the effect of the atom-atom interac- 
tions is no longer negligible and the mean-field nonlinear 
term .giDlV'l 2 m @) starts playing an important role in 
the dynamics; in the following, we shall focus on the ex- 
perimentally most relevant case of a positive scattering 
length a > 0, which means a repulsive interaction among 
atoms in free space. 

In [^6) we discussed the case of a continuous wave (cw) 
beam of coherent atoms incident on a finite lattice; de- 
pending on the detuning of the incident beam with re- 
spect to the frequency of a quasi-bound mode of the lat- 
tice, atom optical limiting or bistability was predicted for 
an incident beam respectively on or above the resonance 
frequency. Here we shall consider the more realistic case 
of a finite atomic wave packet incident on a finite optical 
lattice; its central frequency wq is taken to be close to a 



In the analogy with optics, these quasi-bound modes correspond 
to the resonance peaks of a Fabry-Perot interferometer or, more 
closely, to the cavity modes of a DBR microcavity |2l| . A more 



detailed discussion about quasi-bound modes in optical lattices 
can be found in and in [ll . 
2 Similar oscillationshave been studied in pa] in the case in which 
a large number of closely spaced quasi-bound modes are excited. 
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quasi-bound mode of the lattice and its temporal dura- 
tion much longer than the lifetime of the mode so that the 
frequency spread is narrower than the mode linewidth. 
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FIG. 4: Transmitted pulse shape for Gaussian incident pulses 
on resonance (a) and 0.04 u>bi above (b) resonance with the 
quasi-bound mode of the lattice at u 7 ~ 1.17wBr (solid lines). 
For comparison, same calculations neglecting the interaction 
term (dashed lines) and in the absence of the lattice (dot- 
dashed lines) . Same lattice parameters as in figs.^ and ^| the 
pulse starts at zo/^Br = —1000 with a Gaussian width equal 
to w /£ B r = 480. 
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FIG. 5: Reflected pulse shapes for the same values of the 
physical parameters as in fig.W (black lines). For comparison, 
same calculation neglecting the interaction term (grey lines). 



First we consider the case of an incident wave packet 
with a center-of-mass kinetic energy Hloq exactly on reso- 
nance with a quasi-bound mode (fig^a): in the absence of 
interactions (dashed line), the wave packet is nearly com- 
pletely transmitted without any shape distortion; only a 
small fraction of the pulse is reflected (figj|a). Thanks 
to the resonance condition with the quasi-bound mode, 
the atomic density inside the lattice is significantly larger 
than the one in the incident wave packet. In the pres- 
ence of interaction, the main effect of nonlinearity is to 
blue-shift the quasi-bound mode frequency |11[ and push 



it out of resonance with the incident beam; this negative 
feedback effect, already present in the cw treatment [ fill) , 
gives not only a saturation of the transmission as a func- 
tion of the incident density, but also a significant reshap- 
ing of the wave packet: the higher density part is cut 
and the resulting wave packet is flattened. The negative 
feedback in the transmission obviously corresponds to a 
positive feedback in the reflection. 

In the opposite case of a wave packet incident with 
a kinetic energy higher than the quasi-mode frequency, 
a low density wave packet is nearly completely reflected 
(fig.||b, dashed line). At higher densities, the cw calcula- 
tions in fl6|| have predicted the appearance of atom opti- 
cal bistability effects: several stationary states with dif- 
ferent transmitted intensities can indeed be found for the 
same value of the incident intensity. On such grounds, we 
are able to put forward a physical interpretation of the 
transmitted pulse shape in the presence of interactions 
shown as a solid line in fig.^b: only a small fraction of the 
leading part of the pulse is transmitted since the incident 
frequency is out of resonance with the empty quasi-bound 
mode; moreover, this part of the incident wave packet has 
been accelerated by the repulsive mean-field interactions 
before reaching the lattice and thus pushed even further 
off resonance. For the same reason, the trailing part of 
the pulse has been instead slowed down with respect to 
the central velocity and thus results closer to resonance 
to the quasi-bound mode; when this trailing part of the 
pulse reaches the lattice, the quasi-bound mode starts to 
be effectively populated. The main effect of the inter- 
actions among the atoms in the quasi-bound mode is to 
push its frequency closer to resonance with the incident 
wave packet: the sudden increase in the density of the 
transmitted pulse which is apparent in fig.^b is a direct 
consequence of this positive feedback. This behaviour is 
analog to the jump from the lower to the upper branch 
of the bistability loop which occur in the cw case for inci- 
dent intensities growing beyond the switch-on threshold. 
The shape of the reflected pulse is complementary to the 
one of the transmitted one: in the presence of the non- 
linearity, the reflected pulse shows a dip corresponding 
to the switch-on of the transmission (figjfjb) . 

Provided the interaction energy is much smaller than 
the spacing of the different quasi-bound modes, the trans- 
mission dynamics is mostly determined by a single reso- 
nant mode and the shape of the matter field inside the 
lattice is fixed by the eigenfunction of the mode; this 
guarantees that no modulational instability ^ ^2) can 
take place 3 even in the presence of an effectively attrac- 



3 From a different point of view, this suppression of the modula- 
tional instability can be interpreted in the following terms: if the 
spatial extension of the condensate wave packet is small enough, 
the excitation of the long wavelength modes which are respon- 
sible for the modulational instability can not occur because of 
the finite size of the system; this effect is well known from the 
physics of trapped BECs with attractive interactions, which are 
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tive interaction such as the one which occurs in the case 
of negative mass m c g < valence band atoms for which 
the sign of the effective interaction is reversed with re- 
spect to free space. 

In the absence of spatial confinement, a spatially ex- 
tended wave packet of coherent valence band atoms is in- 
stead subjected to modulational instability: consider for 
example an attractive lattice and an incident wave packet 
with a kinetic energy just below the lower edge of the gap. 
In the present paper we shall limit ourselves to the case 
of a sufficiently weak nonlinearity gm^l 2 ^ Vo in order 
for the band structure of the atomic dispersion not to be 
washed out by the interaction term. Since propagating 
valence band states are available at all spatial positions, 
the pulse is able to penetrate inside the lattice without 
any reflection. As soon as the pulse is in the lattice, 
the modulational instability sets in and the initially uni- 
form envelope starts being modulated with an amplitude 
which grows exponentially in time; at the end, a train 
of short and intense pulses is found (fig.§)). The seed for 
this modulational instability is automatically provided 
by the density modulations which inevitably arise while 
the pulse is entering the non-uniform lattice. In order 
to make the discussion the simplest, this has been as- 
sumed to have an uniform depth in its central part and 
wings at its ends; after the initial preparation phase, the 
pulses therefore propagate through an effectively uniform 
system 4 . 
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stable provided the number of atoms is sufficiently small for the 
effective healing length to be larger than the condensate size . 
4 In the presence of a slowly varying modulation of the lattice 
parameters, no significant coupling of the internal and external 
degrees of freedom of a solitonic pulse is expected to occur if the 
characteristic length of the inhomogeneity is much longer than 
the size of the pulse |24j . 



FIG. 6: Onset of modulational instability: propagation of a 
pulse of valence band atoms (wo/wbf = 0.68, i»o/fer = 240, 
zo/Ibi = —650) across an attractive lattice (Vo/u>bi = —0.4). 
The lattice has a fiat profile in the central \z/£bi\ < 100 region 
and Gaussian (wi/Ibt = 50) wings; its spatial extension is 
indicated by the vertical dashed lines. 

A similar modulational instability is well-known to oc- 
cur in nonlinear Bragg fiber optics Q; a cw laser beam 
traveling in the conduction (valence) photonic band of a 
Bragg lattice with a focusing (defocusing) nonlinear re- 
fractive index is subjected to self-pulsation and therefore 
converted into a train of pulses. Since the nonlinear re- 
fraction index can be interpreted as an effective photon- 
photon interaction |25|, ^6) , the self-pulsing effect is easily 
interpreted as the modulational instability following the 
presence of an effectively attractive interaction. 



V. CONTROLLED BRIGHT GAP SOLITON 
GENERATION 

A regular train of short optical pulses with a well de- 
fined period and intensity has been theoretically shown to 
be produced if the modulational instability of a cw laser 
beam propagating in a lattice is seeded with a weak pe- 
riodic intensity modulation |l(| . The shape of each pulse 
is stable during propagation since the group velocity dis- 
persion is counterbalanced by the optical nonlinearity. 
Following the literature, we shall call these pulses gap 
solitons even if from a strictly mathematical point 
of view they could be classified only as solitary waves, 
since the wave equation ([}]) in a periodic potential is not 
exactly integrable (25) and collisions between two such 
solitons do not exactly preserve the pulse shape ^7j; as 
we shall see in the following, the expression gap soliton is 
however physically justified by the fact that the pulse dis- 
tortion following a collision is generally small. Since the 
first observation of optical gap solitons in 1995, intense 
experimental activity is actually in progress for the gen- 
eration and characterization of gap solitons and Bragg 
modulational instabilities in optical fibers Q. 

In very recent years, solitonic excitations are begin- 
ning to be investigated also in the context of nonlinear 
atom optics [[l8| : dark solitons in the form of stable den- 
sity dips in a otherwise uniform condensate have been 
recently observed |2l|j bright gap solitons Q and modu- 
lational instabilities!™ are actually under intense investi- 
gation; despite the positive atom-atom scattering length, 
such stable atomic pulses can exist inside an optical lat- 
tice thanks to the negative effective mass of valence band 
atoms. In the present section we present a new method 
to generate narrow bright gap solitons starting from a 
wide atomic condensate incident on an optical lattice 

Consider a long coherent matter wave pulse (i.e. a 
Bose-condensed atomic cloud) incident on the same at- 
tractive lattice as in the previous section with a kinetic 
energy just above the lower edge of the gap. At lin- 
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ear regime, such a wave packet is nearly completely re- 
flected so that the interference of the incident and re- 
flected waves creates a standing wave pattern in front of 
the reflection point. Far outside the lattice, the period of 
the pattern is fixed by the wave vector ko of the incident 
wave packet and is therefore of the order of the lattice 
period ts T . Inside the lattice, the standing wave pattern 
originates instead from the interference of a forward and 
a backward Bloch waves: the amplitude of the fast oscil- 
lations follows from the density modulation of each Bloch 
wave function with a periodicity £bi> while the local pe- 
riod of the slower modulation which originates from the 
interference of the two Bloch waves k = &B r ± Afc is equal 
to 2n/Ak. As we approach the reflection point, the Bloch 
waves approach to the band edge, so that Afc — > and 
the period of the modulation is strongly increased with 
respect to the lattice period, although it still remains 
much shorter than the size of incident wave; this effect is 
apparent in the leftmost snapshot of figj^. 
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group velocity v g and the peak density) can be controlled 
by acting on the parameters of the wide incident Bosc- 
condensate. 

Once the pulse has crossed the flat region of the lat- 
tice and has got to its opposite end, the effective mass 
of the atoms becomes positive again and the pulse is im- 
mediately broadened under the combined effect of group 
velocity dispersion and mean field repulsion (sec the two 
last snapshots in fig.0). A proof of the solitonic nature 
of the generated pulse is obtained from a study of its 
collisional dynamics (fig.|^): a pair of such pulses sym- 
metrically generated at the two ends of the lattice collide 
in the middle of the lattice. Their solitonic properties 
result clearly from the fact that their shapes as well as 
the number of atoms contained in each of them are only 
weakly affected by the collision process. The small broad- 
ening of the pulses that can be observed in fig.|| is a sig- 
nature of the fact that the nonlinear wave equation in a 
periodic potential is integrable only in an approximate 
way p7], M. 
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FIG. 7: Controlled bright gap soliton generation from 
a standing matter wave (tJo/wBr = 0.72, wo/lsr ~ 240, 
zo/fer = —650; same lattice parameters as in fig.[]): pulse 
shape snapshots at different times. Vertical dashed lines in- 
dicate the spatial extension of the lattice; the dot-dashed line 
indicate the reflection point at linear regime. 



FIG. 8: Collision process between a pair of gap solitons gen- 
erated at either end of a lattice with the same parameters as 
in figs.^ and propagating along the lattice with opposite ve- 
locities. All the collisional dynamics takes place in the region 
\z/£bt\ < 100 where the lattice profile is flat. 



Provided the nonlinear term is sufficiently small, the 
interactions do not wash out the standing wave inter- 
ference pattern |^9| but simply blue-shift the local band 
edge at the spatial position of the rightmost antinodes; in 
this way, the wave packet frequency is pushed out from 
the local gap and valence band states result available for 
propagation. The resulting short pulses, stabilized by 
the effective attractive interaction, can therefore propa- 
gate along the lattice as solitonic objects. The higher 
the density, the larger the number of pulses for which 
this mechanism is effective and which are then able to 
propagate along all the lattice without being reflected. 
By carefully choosing the density, a single soliton can be 
launched along the lattice (fig-^); unfortunately, given 
the complexity of the nucleation process, it is not physi- 
cally evident how the parameters of the soliton (e.g. the 



Thanks to the effective mass m e ff much smaller than 
the free space mass mo (for the parameters of fig.0 
"ieff = 0.07 mo), the solitonic width is significantly larger 
than the free space healing length at the same value of 
the density; for a typical value of the lattice period of 
the order of 0.5 /xm, the solitonic length turns out of the 
order of 10 /Ltm which is well within the capabilities of 
actual detection systems. The mean-field interaction en- 
ergies required for the observation of solitonic effects are 
in general smaller or of the order of one tenth of the re- 
coil energy WB r : for the most relevant case of 87 Rb atoms 
and 23 Na, this value corresponds to reasonable densities 
of the order of 10 14 cm -3 . The characteristic time for 
the modulational instability and soliton formation pro- 
cesses described in the previous sections is of the order 
of u)B T t ~ 500, which means t ~ 20 ms for 23 Rb atoms 
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and t ~ 3 ms for 87 Na. 

The method here described for the generation of bright 
gap solitons is significantly different from previous pro- 
posals (J] : in our approach the soliton pulse shape origi- 
nates not from the whole BEC cloud, but only from the 
much shorter density bump corresponding to an antin- 
ode of the standing matter wave. This fact allows one 
to obtain short solitons from a wide BEC without the 
need for a dramatic pulse compression under the effect 
of effectively attractive interactions. 



VI. GAP SOLITONS: A SIMPLE ANALYTICAL 
MODEL 

Provided the gap soliton is wide enough, only a narrow 
group of Bloch states around a central wave vector fc so i 
are populated and an accurate description can be ana- 
lytically obtained within the so-called envelope function 
framework: denoting with Uk Bol (z) the Bloch eigenfunc- 
tion at k = fc S oi, we write the wave function ip(z,t) of 
the coherent matter pulse as the product of the slowly 
varying envelope ${z,t) and the fastly oscillating Bloch 
eigenfunction Uk sol (z) 



ij>(z,t) = ip{z,t)u ksol (z); 



(12) 



in the following, the Bloch eigenfunction Uk{z) is assumed 
to be normalized according to 



1 r Br 

— / dz\u k {z)\ 2 



1, 



which corresponds to 

M 2 



\a b \ 2 = l. 



(13) 



(14) 



in the (f , b) basis of (|5|). If its variations of the envelope 
-0(z, t) are slow enough, it can be shown Q that i/>(z, i) 
obeys a simple integrablc nonlinear Schrodinger equation 
(NLSE) 



ih 



dt 



d 2 



2m Q g dz 2 



+ g eff \^(z,t)\ 2 ^(z,t); 



(15) 



in which the effective mass m e ff and the effective coupling 
g e ff depend on the central wave vector k so \; obviously, 
the stability of the solitonic pulse requires m e g < 0, i.e. 
fe so i close to the upper edge of the valence band at /cbi-- 
While an explicit expression for m c g(k) has been already 
given in ([)]), the effective coupling g c s(k) turns out to be 
expressed in terms of the Bloch wave function by 



F(fc) 



1 

fa 



dzgm \u k (z)\ 



(16) 



within the two-mode ansatz (|^), this quantity can be 
rewritten in the simple form 



ges(k) = giu (\df\ 



4|a f | 2 |a b | 2 ) 



(17) 



where a^b are the projections of the Bloch eigenfunction 
on the forward and backward propagating waves. Notice 
that the density modulation of the Bloch wavefunction 
at gap edge (|af| 2 = |ab| 2 = 1/2) makes the effective 
coupling a factor 3/2 larger than the one far from the 
gap, i.e. the free-space one. 



0.06 




FIG. 9: Comparison of the soliton ic p ulse at coBrt = 1700 in 
fig with the analytical prediction ([L8|) for the envelope (solid 
line). Approximate prediction obtained using the band-edge 
values (fegoi = fcer for g e 3 and m e g . 

Under these assumptions, the envelope -0 so i of the soli- 
tonic wave packet has the simple expression pa] 



Aoi(z) = ^ max sech 



with the width £ so i given by 



Z — Vgt 



Csol — ] 

as expected, the size £. 
the healing length £ 



n 2 



m c sg c s \if> 



(18) 



(19) 



,i of the soliton is of the order of 
h/ y/2m c ffg c s iV'maxI 2 . The accu- 



racy of this approximate description is apparent in fig.|^ 
where we compare the numerically obtained wave packet 
with the analytical prediction (18) for the envelope; the 
central wave vector fc so i of the wave packet has been de- 
termined from the group velocity by means of (Q), the 
envelope amplitude i/> max from the peak density of the 
pulse. 



VII. CONCLUSIONS 

In this paper we have theoretically investigated the 
transmission dynamics of coherent matter pulses (such as 
those that can be extracted from Bose-Einstein conden- 
sates) which incide on finite optical lattices; such systems 
are matter wave analogs of the nonlinear Bragg fibers 
currently studied in nonlinear optics. 
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At linear regime (i.e. in the low-density or weak- 
interaction limit), we have characterized the dependence 
of the intensity and shape of the transmitted pulse on the 
velocity and size of the incident condensate in terms of 
the dispersion of matter waves inside the lattice; as in the 
case of light waves in periodic dielectric structures or in 
the case of electrons in crystalline solids, the dispersion of 
matter waves in the periodic potential of optical lattices 
is in fact characterized by allowed bands and forbidden 
gaps. 

The dynamics in the presence of interactions is found 
to be even richer: an interpretation of the numerically 
predicted effects is put forward in terms of familiar con- 
cepts from nonlinear optics, such as optical limiting, op- 
tical bistability, and modulational instability. 

In particular, we have investigated a possible way of 
generating narrow bright gap solitons from a wide in- 
cident Bose-condensate: the modulational instability is 
seeded from the strongly modulated density profile of 
the standing wave which is formed in front of the finite 
optical lattice by the interference of the incident and re- 
flected matter waves. The solitonic nature of the gener- 



ated pulses has been checked from their shape which is 
in excellent agreement with a simple analytical model as 
well as from their dynamical and collisional properties. 

Finally, we have verified that the range of physical pa- 
rameters that is required for the observation of the effects 
predicted in the present paper falls well within the pos- 
sibilities of actual experimental setups. 
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